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A numerical procedure has been developed for the aerodynamic 
force and moment analysis of V/STOL aircraft operating in the 
transition regime between hover and conventional forward flight. 
The procedure specifically treats the interaction between the 
jets and airframe as well as the effect of turbulent flow 
separation on the w:.ng and plain wing trailing edge devices. The 
overall methodology employs three previous] y existing computer 
programs for the calculation of the jet properties and inviscid 
parameters, plus a newly developed method for predicting wing 
viscous effects. 

The trajectories, cross sectional area variations, and tiass 
entrainment rates of the jets are calculated by the Adler-Baron 
Jet-in-Crossflow Program. The inviscid effects of the 
interaction between the jets and airframe on the aerodynamic 
properties are determined by use of the MCAIR 3-D Subsonic 
Potential Flow Program, a surface panel method. In addition, the 
MCAIR 3-D Geometry Influence Coefficient Program is used to 
calculate a matrix of partial derivatives that represent the rate 
of change of the inviscid aerodynamic properties with respect to 
arbitrary changes in the effective wing shape. 

For each baseline configuration, the calculated quantities 
from the second and third programs above establish an input file 
to the MCAIR Stalled Wing Analysis Program (SWAP). The purpose 
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of SWAP is to calculate the surface pressure distribution, forces 
and moments on the aircraft in the presence of wing viscous 
effects including turbulent flow separation on the wing and 
unslotted wing trailing edge devices. The viscous-inviscid 
interaction is explicitly modelled by a first order coupling 
between mathematical expansions of the inviscid and viscous flow 
methods. SWAP can be employed for wing alone geometries, 
wing-fuselage combinations, as well as wing-fuselage- jet 
configurations . 

The complete calculation procedure is described with mathema- 
tical formulations presented for the Stalled Wing Analysis Pro- 
gram. Example solutions are presented that demonstrate the 
accuracy, numerical stability, and efficiency of the jet-aircraf t- 
viscous interaction methodology. 
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INTRODUCTION 

An automated method is presented for predicting the aero- 
dynamic performance of complete V/STOL aircraft operating out of 
ground effect in the transition regime between hover and conven- 
tional forward flight. In transition, the large jet injection 
angles, low aircraft velocity, relatively low Reynolds number, 
and high trailing edge flap deflection .1 lead to strong viscous 
aerodynamic-propulsion interactions that have a substantial 
effect on aircraft forces and moments (Figure 1). 



a deg Cm AERo'° 103 C ' 

FIGURE i. V/STOL TRANSITION: EFFECT CF JETS ON AERODYNAMIC FORCE AND MOMENT 
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f The specific flow phenomena that produce these large inter- 

actions are presented in Figure 2. Although the blockage and 
mass entrainment of the jets have a substantial effect on the 
surrounding aerodynamic flow field, the effect of the airframe on 
the properties of the jets-in-crossflow is typically very weak. 
This fact allows the jet centerline trajectory, cross-section 
geometry, and entrainment rate to be calculated by an independent 
(decoupled) jet-in-crossflow method. 


SEPARATED WAKE ON CONTROL SURFACES/WING 



JET BLOCKAGE 
AND ENTRAINMENT 


FIGURE 2. PROBLEM DESCRIPTION: V/STOL TRANSITION FLOW FIELD 
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If the surrounding aerodynamic flow field were incompress- 
ible and inviscid, a surface panel method could be applied to 
calculate the induced forces and moments on the airframe surface 
(Reference 2). Panelled jets with surface suction would 
represent the effects of jet bLockage and entrainment. However, 
in the transition mode there can be large regions of separated 
flow associated with highly deflected control surfaces or high 
angles-of-attack . Panel methods are appropriate for analyzing 
the flow in the presence of regions of separation if the panelled 
geometry includes the viscous displacement surface. The problem 
is to solve for the displaced surface; the difficulty is that the 
solution is dependent upon strong viscous-inviscid interactions. 

In the past decade several methods have been developed for 
predicting the incompressible aerodynamic forces and moments of 
two dimensional airfoils in the presence of turbulent flow 
separatio (References 3-7). Most of these methods empirically 
model the turbulent flow separation region by constant pressure 
from the separation point to the airfoil trailing edge with some 
method for closing off the wake downstream. This bypasses the 
need to perform complex computations within the viscous separated 
flow region. An inverse two-dimens lonal potential flow method is 
used to "design" the streamline contour that satisfies the 
constant pressure conditions on the separated wake. 

With the development of three-dimensional wing inverse poten- 
tial flow (or design) methods ^References 6-11), the expectation 
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g was that finite winy viscous-inviscid interactions includiny 

turbulent flow separation would be f orthcominy . However, 

reliable methods had not been developed for several reasons. 
First, most of the two dimensional viscous-inviscid interaction 
procedures utilize simple iteration between distinct viscous-only 
and inviscid-only calculation steps. This procedure can be 
numerically unstable and often requires man- in- the- loop 
corrections to insure accuracy and numerical stability. Second, 
accurate three-dimensional conventional panel methods and design 
methods are expensive to use. Iteratiny between a panel method, 
a viscous calculation method, and a design method for five or 
more iterations for each anyle-of-attack can be prohibitively 
expensive. Third, until recently the design methods suffered 
from numerical instabilities that precluded the accurate and 
efficient coupling to viscous methods. All of these difficulties 
should be dealt with simultaneously when considering the finite 
wing viscous-inviscid interaction problem. 

Three recent advances have overcome the major difficulties 
with three dimensional viscous-inviscid interaction predictions. 
In 1980 the Mct'kinnell Aircraft Company IMCAIR) Stalled Airfoil 
Analysis Program (Reference 7) was developed for the reliable 
prediction of 2-D airfoil forces ana moments in fully attached or 
turbulent separated flow In 1981, the MCAIR 3-D Perturbation 
Analysis Method (Reference 12) was developed under contract to 
NASA, Langley Research Center, for the accurate and inexpensive 
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inviscid solution corresponding to arbitrary small perturbations 
to wing-fuselage geometries. In 1982, the MCAIR 3-D Subsonic 
Wing Design Program (Reference 11) was developed under contract 
to NASA, Langley Research Center, for the accurate, inexpensive, 
numerically stable design of wing-on-fuselage geometries 
corresponding to prescribed pressure distributions . 

This report presents a method developed at MCAIR for the 
analysis of jet-airframe wing viscous-inviscid interactions 
including turbulent flow separation on the wing or plain win r j 
trailing edge devices. The method, the Jet-Aircraft-Viscous 
Interaction (JAVI) Program, represents the synthesis of the three 
advances that make possible accurate, numerically stable 
solutions to the 3-D viscous flow problem. Furthermore, by 
incorporating complex aircraft configurations through the use of 
the panel method, the approach is applicable to wings, 
wing-fuselage combinations, and wing-fuselage- jet configurations. 

This report presents the complete calculation procedure for 
treating jet, airframe. and wing viscous interactions. The four 
different computer programs employed are discussed. Mathematical 
formulations for the wing viscous-inviscid interaction analysis 
are presented with example solutions that demonstrate the 
accuracy and stability of the method. Conclusions and 
recommendations are presented. 
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CALCULATION PROCEDURE FOR JET-AIRFRAME -VISCOUS INTERACTIONS 


The computation of the aerodynamic performance of V/STOL air- 
craft in transition out of ground effect requires the use of four 
computer programs: 

la) Adler-Baron Jet-In-Crossflow Program (JICP; Ref 13) 

(b) MCAIR 3-D Subsonic Potential Flow Program (SPFP; Ref 
12 ) 

(c) MCAIR 3-D Geometry Influence Coefficient Program (GICP; 
Ref 12) 

(d) MCAIR 3-D Stalled Wing Analysis Program (SWAP) 

It is assumed the the flow is incompressible and that the follow- 
ing quantities are known: (1) airframe geometry, (2) the free- 

stream angle of attack and unit Reynolds number, (3) the engine 
inlet ma; s flow, and (4) uniform nozzle exit conditions 
(injection angle and jet to freestream velocity ratio). 
Furthermore, it is assumed that the jets are incompressible and 
turbulent, and any trailing edge device is unslotted. 

The calculation procedure is depicted in Figure 3. The 
first step is to determine the jet properties under the assump- 
tion that airframe influence on the jets is negligible. Each jet 

Wj 

injection angle, *j, and jet-to- freestream velocity ratio, r- -» 
are converted by JICP to jet centerline trajectory, cross sec- 
tional geometry, and mass entrainment rate. The jet geometry is 


6 



ORIGINAL PAGE JS 

OF POOR QUALITY 


panelled and connected to a panelled representation of the 
baseline aircraft geometry (aircraft without viscous displacement 
thickness). The mass entrainment of the jet is represented by 
suction distributed over the jet panels. 



JET GEOMETRY AND 
MASS ENTRAINMENT RATE 



(d$|/dZ|) a m q. 90 * 


VISCOUS GEOMETRY. 

a. 


3-D STALLED SWEPT 
WING PROGRAM 
(ITERATION REQUIRED) 

1 1 



FIGURE 3. CONFIGURATION VISCOU® NVISCIO ANALYSIS PROCEDURE 
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The baseline configuration is converted to a distribution of 
inviscid surface pressure and an aerodynamic influence 
coefficient matrix by SPFP. The output from this program is used 
as input to GICP, which calculates a matrix of partial deriva- 
tives. Each element of this matrix represents the rate of change 
of the inviscid perturbation potential on the surface with 
respect to arbitrary changes in the wing displacement. GICP is 
executed twice corresponding to angles-of-attack of zero and 
ninety degrees. 

The output from SPFP and GICP are used as input files to 
SWAP, which usee the inviscid perturbation matrices from GICP to 
extrapolate for the pressure solution corresponding to the 
desired angle-of-attack and calculated wing displacement surface. 
The wing effective shape is determined iteratively along with the 
distribution of pressure over the aircraft surface. Forces and 
moments are determined by integration of pressure over the 
aircraft and surface shear stress over the wing. 

The remainder of this section is devoted to a description of 
the four programs. More detailed descriptions of the first three 
programs are available in the literature (References 12 and 13), 
and the mathematical formulation for SWAP is presented in the 
next section. 
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The Adler-Baron Jet-In-Crossflow Program (Reference 13) con- 
verts the jet injection angle and jet-to-f reestream velocity 
ratio to the jet centerline trajectory, cross-sectional area 
variation, and mass entrainment rate. The method is applicable 
to incompressible turbulent jets and is based on two integral 
momentum equations. One of these is in the direction 

perpendicular to the jet centerline, and the other is in the 
direction parallel to the jet centerline. The mass entrainment 
rate, needed to solve these equations, is based on a linear 
combination of straight jet entrainment (no crossflow) and vortex 
pair entrainment. Integration of the two momentum equations 

along the jet centerline requires a knowledge of the velocity 
distribution over the jet cross section. The contour of jet 

velocity at the boundary between the jet and freestream is 

determined by calculating the displacement of two-dimensional vor- 
tices initially seeded on the nozzle exit. The vortices simulate 
the shear layer between the jet and surrounding flow. The inner 
contours of progressively higher velocity are generated by 

solving Poisson's equation, which is an empirical mathematical 
model. A detailed discussion of the procedure is found in 

Reference 13. 

The Jet-In-Crossflow Program is restricted to one isolated 
jet issuing from a flat plate into freestream. Due to the weak 
effect of the airframe on the jet, it is reasonable to calculate 
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the jet properties in isolation. For multiple jets, the program 
is executed once for each jet. Furthermore, for tandem jets, as 
is the case for the YAV-8B , the upstream jet exerts a large 
influence on the downstream one, and the jets tend to coalesce. 
For multiple jets with interaction and coalescence, the method of 
Wooler (Reference 14) is used to determine the blockade effects 
of the upstream jet on the downstream jet. The merged single jet 
properties are determined to a first approximation by simply 
combining the effects of the individual jets without coalescence. 


A demonstration of the accuracy of JICP is shown in Figures 
4-6, where the predicted jet centerline trajectory, cross- 
sectional area ratio, and entrained mass flux are compared with 
experiments (Reference lb). The agreement is good for the cases 
examined. 


Z' 



FIGURE 4. LOCUS OF PEAK JET VELOCITY CENTERLINE TRAJECTORY 

90* INJECTION ANGLE 
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FIGURE 5. JET CROSS-SECTIONAL AREA RATIO COMPARISON 
90* INJECTION ANGLE 



FIGURE 8. ENTRAINED MASS FLUX COMPARISON 

90* INJECTION ANGLE 
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Subsonic Potential Flow Proyram 

The MCAIR 3-D Subsonic Potential Flow Program (Reference 12) 
is a surface panel method that is based on the combined source- 
doublet distribution of the classical third identity of Green 
(Reference lb). The advantages of this combined source-doublet 
distribution for a surface panel method are well documented 
(Reference 17). The mathematical formulation employs a constant 
source distribution and a quadratic doublet distribution on each 
panel, where the solution singularity strengths are determined by 
satisfying indirect internal perturbation potential boundary 
conditions (Reference 18). The flow velocity on each panel is 
then established from local velocity-singularity relationships 
associated with Green's third identity, instead of direct 
summation of the influences of the singularities on all the 
panels. The method is formulated with a G^thert type 
compressibility correction; however, this option is not used in 
the present method. A detailed discussion of the method is found 
in Reference 12. 

This panel method is used to calculate the airframe inviscid 
pressure distribution induced by the aircraft forward velocity 
plus the blockage and mass entrainment effects of the jets. A 
demonstration of the accuracy of the panel method is presented in 
Figure 7, where the YAV-8B wing pressure distribution in conven- 
tional flight (.30 , o.4* angle-of-attack ) is compared to wind 

tunnel data (Reference 1). The ability of the panel method to 
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model geometric details is reflected in the accuracy of the lower 
surface pressures at 25% semispan, where the effect of the 
protruding nozzles is substantial. Wind tunnel surface pressure 
data for the YAV-8B in V/STOL transition is not available for 
comparison with analytical predictions; however, detailed 
pressure data are available for a single circular jet emanating 
from a flat plate into a uniform crossflow (Reference 19). In 
Figure 8, the geometry for the panelled jet was calculated by the 
Adler-Baron method. In the surface panel method, the effect of 
the mass entrainment is simulated by distributed normal 
velocities on the jet panel surfaces. The pressure distribution 
on the plate calculated by the panel method agrees well with 
experimental data, except in the wake of the jet (Figure 9). 
Although the separated wake of a jet-in-crossflow represents an 
important aerodynamic problem that requires further 
investigation, it is beyond the scope of the present effort. The 
good agreement for this jet-plate combination demonstrates that 
the treatment of the jet decoupled from plate is satisfactory. 

Geometry Influence Coefficient Program 

The formulation for this program is similar to the conven- 
tional panel method with one major exception. In each step of 
the solution process, the derivative of each quantity with 
respect to geometric coordinates is established, rather than the 
quantity itself. This process culminates in the generation of a 
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matrix of partial derivatives of potential with respect to arbi- 
trary geometric perturbations. The complete mathematical 
formulation is described in Reference 12. 


QUALITY 


c B 



p«n#t method 

A Q Wind tunnel data (Sat 1) MS PANELS (TAIL-OFF) 
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FIGURE 7. YAV-8B AT 6.4* ANGLE-OF-ATTACK >N CRUISE MOOE 

Mqo =0.50 


The purpose of establishing the inviscid perturbation matrix 
is to eliminate the repetition of computationally expensive steps 
corresponding to a series of small arbitrary geometry perturba- 
tions to a given baseline configuration. By running this program 
twice at angles-of-attack of zero and ninety degrees and using an 
automated extrapolation procedure incorporated in SWAP, the dis- 
tribution of surface pressure, corresponding to arbitrary angle- 
of-attack and small geometry perturbations to the baseline, can 
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be calculated at small expense. The method ia very accurate for 
large changes in wing camber and twist, characteristic of wing 
viscous displacement effects. This extrapolation procedure is 
called the Perturbation Analysis Method (Reference 12 ) . 
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As illustrated in Figures 10-15, the Perturbation Analysis 
Method is accurate for surprisingly large geometry perturbations. 
The low aspect ratio wing of Figure 10 with constant NACA 0012 
section geometry was selected as a baseline panelled 
configuration. The section geometry was then perturbed twice, 
first to form a supercritical wing and second to form a fiyhter 
wing (Figure 11). Computed results for the perturbed geometries 
(Fiyures 12-15) are nearly identical to the virtually exact 
solutions obtained using a conventional surface panel method, but 
the computing expense is more than an order of magnitude less . 

Stalled Wing Analysis Program 

The objective of the Stalled Wing Analysis Program is to 
determine the wing effective shape (viscous displacement surface) 
and the resulting distribution of surface pressure over the air- 
craft. JICP and SPFP are each executed once for each case to be 
analyzed. GICP is executed twice. The data needed to run SWAP 
are stored as permanent input files. The analysis of the 
geometry over a range of anyles-of-attack and different Reynolds 
numbers may be then accomplished by SWAP without having to rerun 
any of the first three computer programs. 
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FIGURE 10. BASELINE WING PANELING FOR DEMONSTRATION 
OF PERTURBATION ANALYSIS METHOD 
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FIGURE 11. TEST CASES FOR THE PERTURBATION ANALYSIS METHOD 


The procedure in SWAP is to (1) extrapolate the inviscid 
pressure distribution for the anyle-of-attack and latest estimate 
for the wing effective shape, (2) calculate the viscous flow 
parameters and the difference between the denired conditions that 
define the wing effective shape and the actual calculations, (3) 
linearly couple the inviscid and viscous parameters to determine 
the changes in the effective shape necessary to simultaneously 


satisfy the desired 

conditions to 

first 

order. 

(4) update 

the 

effective shape, and (5) repeat 

steps 

(1) 

through (4) 

to 

convergence . The 

forces and 

moments 

are 

determined 

by 


integration of pressure over the aircraft surface and wall shear 
stress over the wing. The mathematical formulation and example 
solutions are presented in the next two sections. 
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FIGURE 12. PRESSURE DISTRIBUTION OF WING A AT 5° ANGLE OF- ATTACK 

SUPERCRITICAL 
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0 0.2 0.4 0.6 0.8 1.0 

FIGURE 15. FORCE AND MOMENT DISTRIBUTION OF WING B AT 5* ANGLE OF ATTACK 

FIGHTER 
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SWAP is operational on the MDC Cyber 17 5 and is accurate and 
reliable for fully attached flow. The method requires further 
development for separated flows as expected reliability has not 
been achieved. Prediction at angles-of-attack below maximum lift 
are satisfactory; however, as maximum lift anyle-of-attack. is 
approached, the method has difficulties with both convergence and 
accuracy. These problems are being investigated with solutions 
expected in the near future. 


MATHEMATICAL FORMULATION FOR STALLED WING ANALYSIS PROGRAM 

The Stalled Wing Analysis Program (SWAP) is formulated to 
allow for turbulent separation on the wing and one full span 
plain wing trailing edge device. Short bubble type laminar 
separation is allowed and is treated as the point of transition 
to turbulence. Long bubble laminar separation or bubble burst is 
not treated. 

The method is based on the hypothesis that there exists a 
viscous displacement surface (wing effective shape) on which the 
pressures calculated by potential theory agree with the viscous 

pressures on the surface of the wing. The formulation requires 
(1) a set of conditions, called the theoretical model, that 
specify this displacement surface, (2) numerical methods for 
determining the surface pressure and boundary layer growth, (3) a 
specialized, discretized representation of the geometry and (4) a 
calculation procedure for determining the geometry that satisfies 
the theoretical model. Each of these topics is discussed. 
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Theoretical Model for Wing Effective Shape 


The set of conditions that specify the winy effective shape 
is depicted in Figure 16. The viscous displacement of the 

laminar boundary layer and the forward half of the length of the 
turbulent boundary layer on both upper and lower surfaces is 
neglected. The transition of the boundary layer from laminar to 
turbulent flow is considered instantaneous and is determined from 
one of three empirical correlations for mainflow transition, 
crossflow transition, and laminar separation. For the aft half 
of turbulent attached flow, the effective wing shape is the 
conventional representation of the wing plus boundary layer 
displacement thickness, 5*. The location of the turbulent 
separation line is based on a critical value of the mainflow boun- 
dary layer shape factor, H, which is representative of vanishing 
wall shear stress in the mainflow direction. 

The separated viscous wake model is based on empirical obser- 
vations. First, the pressure within the separation region along 
the wing surface in the chordwise direction is approximately 
constant. Downstream of the wing trailing edge, the viscous wake 
is divided into two regions, referred to as the fore and aft 

trailing viscous wakes. The fore wake is allowed to have 

thickness and camber, whereas the aft wake is very thin having 
only camber. The wake camber line is specified by the condition 
that the wake cannot support a force, ^\C p = U. The thickness 

distribution in the fore wake is determined by requiring the 
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/" chordwise pressure recovery on this surface to be linear from the 

value at the wing trailing edge to approximately freestream at 
the end of the fore wake. 




FIGURE 16. THEORETICAL MODEL FOR VISCOUS-INVISCID INTERACTIONS 


v 
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The inviscid trailing wake is modelled as a sheet of 
constant streamwise vorticity extending far downstream. The 
effect of the roll up of the flap ind wing tip vortices is 

considered to be small compared to the effect of viscous 

separation and is not modelled. 

Selected Numerical Methods 

The numerical methods were selected based on accuracy and 
numerical stability. As mentioned in the previous section, the 
pressure solution is determined by incorporating the logic of the 
Perturbation Analysis Method (Reference 12). In addition to the 
accuracy and efficiency previously discussed, this method offered 
the advantage of providing a mathematical formulation for 

calculating a matrix, whose elements correspond to the rate of 
change in the velocity at each panel center with respect to 

arbitrary displacements to the panel cornerpoints . This matrix 
is necessary to linearly couple the inviscid and viscous 
conditions of the theoretical model. 

The selected methods for determining the boundary layer 
growth on the wing surface are depicted in Figure 17. Along the 
attachment line, the laminar flow, transition, and turbulent flow 
are predicted by the methods due to Rosenhead, Pfenninger, and 
Smith (References 20, 21, and 22), respectively. The laminar 

flow characteristics are predicted by the Cooke method (Reference 
23). Mainflow and crossflow transition are determined by 
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empirical correlations due to Michel (Reference 24) and Gross 
1 Reference 25;. The turbulent flow characteristics are predicted 
by the method due to Smith (Reference 2b). The separation loca- 
tion is determined by the mainflow shape factor exceeding a value 
of 2.0. 



FIGURE 17. SELECTED BOUNDARY LAYER METHODS 


All of these methods are integral boundary layer methods, 
or, in the case of the transition prediction methods, analytical 
representations of empirical curve fits. The empirical curve 
fits were selected primarily because efficient, accurate, analyti- 
cal methods are not available for predicting three-dimensional 
transition. The integral boundary layer methods were selected 
for their accuracy and numerical stability. While not as accu- 
rate as finite difference methods in predicting detailea boundary 
layer parameters, the integral methods are entirely sufficien - 
for calculating tne global effects of viscous-inviscid 
interactions on the total forces and moments. 


\ 
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Geometric Representation of Wing Effective Shape 


A right-handed Cartesian coordinate system is employed by 
the inviscid calculation methods. X is taken in the chordwise 
direction from leading to trailing edge with positive Y pointing 
outward from wing root for a right handed wing. All geometry 
indexing for the wing proceeds from the viscous wake trailing 
edge along the upper surface to the lower surface viscous wake 
trailing edge. Changes due to viscosity to the baseline 

configuration ( jets-fuselage-wing without viscous effects) are 
allowed only on the wing surface. Thus, even though the presence 
of fuselage and jets is accounted for by the panel method 
calculation, the viscous-inviscid analysis procedure treats the 
wing as if it were isola^d. 

The geometric representation of the wing effective shape 


must 

'.nterface 

with 

the 

Perturbation 

Analysis 

Method (Reference 

12), 

which is 

used 

to 

compute the 

inviscid 

pressures. This 


method represents the effective shape as displacements in the Z 
direction from the baseline (Figure 18). These displacements are 
applied at the geometric defining stations, the panel corner 
points. The wing geometry is represented by panels with lines 
connecting corner points in the chordwise direction called 
defining span stations. However, all span stations cannot be 
independent since the number of prescribed pressures in the 
theoretical model must exceed the number of geometric 

perturbations for numerical stability (Reference 11). The even 
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numbered span stations on the wing are considered linearly 
dependent upon the two surrounding odd numbered span stations . 



FIGURE 18. GEOMETRIC DESCRIPTION OF WING EFFECTIVE SHAPE 


The viscous flow computations employ a body iitted non- 
orthogonal coordinate system (Figure 19), where X' is the surface 
distance from the upper surface viscous wake trailing edye, and 
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w* 



FIGURE 19. VISCOUS FLOW COORDINATE SYSTEM 


Y* is parallel to the leading ed^e sweep. The wirvg surface is 
considered flat with 2' taken in the positive Z direction on the 
upper surface and in the negative Z direction on the lower sur- 
face. The surface distance to the panel cornerpoints , X', is 

ijiven by: 


l-l 

X*i - Z ClXj+1 - X : ) 2 + IZ-j+i - Z-,) 2 ] U) 

j-n 
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where Ji, is the first panel cornerpoint on thi* span station. il* 
is determined by first dividing each of the independent span sta- 
tions into fourteen possible viscous rations l Figure 20) 
depending on the physics of the flow. Displacement of the win* 
effective shape is neglected in regions o-9. The displacement of 
the remaining regions is represented by cubic polynomials based 
upon the values IX*. I', t* 1 ) at the region endpoints iFi>jure .1), 

J' • • 

where ' is measured relative to baseline 1 t* ‘ - tan"‘ j~-r). °hce 
the latest set of independent region endpoint conditions has been 



REGION NUMBER 

VISCOUS REGION 

1 

AFT UPPER SURFACE VISCOUS TRAILING WAKE 

2 

FORE UPPER SURFACE VISCOUS TRAILING WAKE 

3* 

AFT UPPER SURFACE SEPARATED REGION 

4** 

FORE UPPER SURFACE SEPARATED REGION 

5 

AFT UPPER SURFACE TURBULENT FLOW 

fl 

FORE UPPER SURFACE TURBULENT FLOW 

T 

UPPER SURFACE LAMINAR FLOW 

A 

LOWER SURFACE LAMINAR FLOW 

9 

FORE LOWER SURFACE TURBULENT FLOW 

10 

AFT LOWER SURFACE TURBULENT FLOW 

If* 

FORE LOWER SURFACE SEPARATED REGION 

12* 

AFT LOWER SURFACE SEPARATED REGION 

13 

FORE LOWER SURFACE VISCOUS TRAILING WAKE 

14 

AFT LOWER SURFACE VISCOUS TRAILING WAKE 


*£%«•!• only i» *#p*f#tton •» Ah# of 'Up hmg# no# 

only if ««p#i # fton • 01 ••#nt 


FIGURE 20. VISCOUS REGIONS ALONG SPAN STATIONS 
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j 


determined by the viscous-inviscid interaction procedure, the 
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a l' a 2' a 3' a 4 " functions of ( v ‘ A . 2 ‘ A , 0‘ A , X * B , Z ' B , 0' B ) (4) 

The baseline configuration is allowed to have surface slope 
discontinuities at the flap hinge line, if flaps are present, and 
at the wing trailing edgo. However the wing effective shape 
must be continuous at all poiits except where the fore wake joins 
the aft wake (Figure 22 ) The baseline discontinuities are 
converted to a continuous displacement surface by requiring: 


where 


e 'a ■ e‘ * e'D 

i?' ■ angle allowed to change during perturbation 

0’ D » discontinuous angle in baseline 


Thus 0' D accounts for the kink in the baseline shape. 

At the fore wake trailing edge the effective shape is 
allowed to be discontinuous (See Figure 22). The downstream 
angle, 0^, is required to be the bisector of the two upstream 
angles, 0’wu and 0'wL,: 


WM 


4 ’wu 


+ 

2 


WL 


lu) 
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The discontinuity is artifically introduced by requiring 9' and 
9 'd in equation (5) to be (on the upper surface): 

9 ' * 9'wm (7) 

9 ‘d * 9 'WU “ 9 ' WM (8) 


Similar equations are employed for the lower surface. 

On the basis of numerical studies, it is reasonable to 
represent the length of the fore wake as a linear function of the 
thickness of the displacement surface at the wing trailing edge 
(See Figure 23). The curve shown in Figure 23 was generated from 
numerous computer solutions for different airfoils and Reynolds 
number by the Stalled Airfoil Analysis Program (Reference 7). 
This program solves for this length by requiring the velocity to 
be linear from the airfoil trailing edge value to approximately 
freestream at the wake trailing edge. However, in the present 
method, this procedure cannot be used due to the surface slope 
discontinuity in the effective shape at this point. 

The turbulent boundary layer is divided into two regions, 
which are of equal length if flaps are not deflected or 
separation occurs upstream of the hinge line. Otherwise, the 
juncture between the two regions is fixed at the flap hinge line. 
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FIGURE 22. BASELINE SHAPE SHOWING SLOPE DISCONTINUITIES 



DISPLACEMENT THICKNESS AT WING TRAILING EDGE (Z'/C) 

FIGURE 23. LENGTH OF THICK VISCOUS WAKE AS A FUNCTION OF DISPLACEMENT 
THICKNESS AT WING TRAILING EDGE (2 D) 
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The polynomial representation of the displacement surface in 
separated regions is made quartic in nature by the addition of 
another unknown displacement Z'^, at the midpoint of the separa- 
tion region. For these regions the displacement at panel 
cornerpoints is given by: 


iZi = AZ Ci + lb(Z’ M - Z* MC )(5i 4 - 25i 3 +5i 2 ) (9) 

where AZ^ is given by equation (2) and Zm^ is given by equation 

(2) for » .5. Only one separation region is normally used on 

either the upper (P*4) or lower surface (P*ll). However, when 

separation occurs ahead of the flap, both separation regions are 

employed with the endpoint between the two fixed at the hinge 
line . 


Incorporation of the above assumptions and approximations 
and noting that the length of the trailing viscous wake is fixed, 
results in a relatively small set of independent unknowns, 
depicted in Figure 24 for the case of separation aft of the flap 
hinge line on the upper surface. The number of unknowns per span 
station can be as low as ten for no separation without laminar 
flow and can be as high as twenty for separation ahead of flap 
hinge line with laminar flow. 
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FIGURE 24. INDEPENDENT UNKNOWNS FOR A SPAN STATION WITH UPPER SURFACE 
SEPARATION AFT OF FLAP HINGE LINE 


The set of independent unknowns across the wing can vary 
from one independent span station to another. The complete set 
of unknowns for the wing is stored in an array, called the 
"g-array " : 


— — 


— - 



91 

92 


z ‘ WE 
0 'WE 

0 ’WL 

► 

Independent 
Span * 1 



# - 

* 

2 





NSPANI 

*N un k 


- 




The determination of the g-array that results in most closely 
satisfying the conditions of the theoretical model uniquely 
specifies the desired effective shape. 
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The calculation procedure employed in the Stalled Winy Analy- 


sis Program is depicted in Figure 25. It is assumed that the geo- 
metric configuration to be analyzed has been previously analyzed 


by the MCAIR 3-D Subsonic Potential Flow Program and MCAIR 3-D 


Geometry Influence Coefficient Program. Thus, two matrices, 


corresponding 

to zero 

and 

ninety 

degree angles' 

-of-attack, are 

available 

as 

input 

files . 

The 

elements 

of 

these matrices 

represent 

the 

first 

order 

rate 

of change 

of 

potential with 


respect to arbitrary changes in the AZ's on the independent span 
stations on the wing. 

The user provides as input the freestream Reynolds number 
based on reference chord, velocity magnitude, and angle-of-attack 
range to be analyzed. In addition, the viscous geometry is 
chosen by specifying the number of calculation points and 
distribution for each viscous region. An initial guess for the 
location of transition is chosen. The user also has options to 
fix transition at the input value, specify that the boundary 
layer flow on the wing be fully turbulent, or fix separation at 
an input location. From these input values, the g-array is 
initialized. 
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FIGURE 25. CALCULATION PROCEDURE OF STALLED WING ANALYSIS PROGRAM 
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» Equations (2) - 19) are employed to determine the set of Az^ 

from the g*. The magnitude of the velocity at each panel center 
i is determined from Green's third identity (Reference lo ) after 
the potential has been extrapolated from: 


■ cosa( $0^ 


NKS 
+ Z 
3-1 


, 


*i , a- 0 . AZ 3 ) * = (TzT^.9u* AZ 3 


NKS 

4 

3" 




( 11 ) 


where and are the baseline potentials corresponding co 

a »■ U* and 90* from SPFP and NKS is the number of independent 
perturbations, AZ. The local velocity direction is determined 
from panel orientation, which allows the determination of the 
velocity components in the X, Y, and Z directions (W x , Wy, W z )i. 
The tangential components of the inviscid velocity, U and V in 
the chordwise and spanwise direction are determined from: 


Ui 

v i 


W x x N Zl - W z x N x x 
(N Xi 2 ♦ N 2 i 2 )l ^ 5 


( 12 ) 

(13) 


where (Nx, Ny. Nj ) i are the components of the unit normal vector 
for the ith panel. This definition of U provides a negative 
value of U on the lower surface of the wing aft of the attachment 
line. Thus, the location of the attachment line, the boundary 
between the upper and lower surface laminar regions, can be deter- 
mined from the point at which U is equal to zero along each span. 
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The pressure coefficient at the panel center, Cp x , is calcu- 
lated from Bernoulli's equation for steady flow: 


c Pi - 1 - 


u i 2 ♦ v i 2 
W 2 


(14) 


In addition, the rates of changes to U and V due to arbitrary 

changes in the Zj can be calculated from equations (11) - (13) 

and the velocity direction equation. These coefficients are 
dUi dV| 

stored as and ^ . These matrices are converted to rates of 

change with respect to changes in the g-array by: 


dUi 

dUi 

dZj 

dik 

dZ^ 

* dg* 

dVi 

. <™i 

dZ j 


dZj 

dgk 


(15) 


vlt>) 


dZ- 


The is determined by differentiation of equations (2), (3) 


(4), and (9). 


The next step in the calculation procedure is to determine 
the coiditions at the viscous calculation points, which are taken 
along the panel centers. The tangential velocity components, IK 
and V., at these points are determined by cubic interpolation in 
the X' direction on Uj_ and Vj,. Thus, 

Ui-i * b 2 U; ♦ I 3 Uj^i ♦ b 4 Ui^-2 * ^5 ^'vi 117) 


\ 
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where 

t>l,b 2 .b 3 #b 4 ,b 5 ■ coefficients based on interpolation routine 

• functions of X' v? (lb) 

and the location of the matching point is between Ui ana 

^i+l* From equations (17) and (18)i 



&1 


du i-i , , dU i , , du i + i 

-35T * b2 55J + -d5T 


dU 


b 4 


i+2 


d 9k 


3 b)c 

* (b6 + k L «'»<' 


dX' V; 

d 9k 


(19) 


Similar calculations are performed for V. 


S 

V 


The values of X' v and Z' v are determined from the y-array: 


X 1 


v l 


Ay(lT+ y(2) * 'i (x 'b l " x 'a l ) 
Ay(i)+ y( 2 ) ^ ’ (x b r x ’ a r^ 


(20) 


z* 




Ay(l)+Ay(2) 




z' . ♦ 

L Ay(l)+iy(2) 


z ' 

- R (21) 


where Ay(l), and Ay(2) are defined in Figure 2b, L and R refer to 
left and riyht independent span stations surroundiny the panel 
centers being calculated, (, is the input value of £: 




X' - X‘ A 

— i 

X*B - X* A 


V 


(22) 
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A and B refer to region endpoint conditions, 
determined from equations (2), (3), (4) and (9). 



are determined by differentiating equations (20) and (21). 



FIGURE 26. SCHEMATIC OF BOUNDARY LAYER CALCULATION POINTS 


Once the inviscid condition* are established at the boundary 
layer calculation points, the viscous parameters can be deter- 
mined. In addition to the inviscid parameters previously 
discussed, the boundary layer calculation routines *. -“quire grad- 
ients in the X and Y directions of the velocity components U and 
V: U x , Uy, V x , Vy. The X direction gradients are established by 
differentiating equation (17) for L’ x and the analog for V x . The 
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/ Y direction gradients are calculated from differentiation of a 

quadratic curve fit in the Y direction. These terms are also dif- 
ferentiated with respect to g to establish the perturbation 
matrices . 

The boundary layer calculation methods are based on the hypo- 
thesis that the viscous flow can be divided into components in 
the mainflow and crossflow direction with the components behaving 
much like two-dimensional flow. A sketch of the boundary layer 

coordinate system is depicted in Figure 27. The angle between 
the projection of the external streamline onto the surface and 

the X' axis, denoted Y, is the direction of the mainflow. The 
crossflow direction is taken perpendicular to this direction in 
the positive Y‘ direction. The angle between the mainflow 
direction and the limiting wall streamline, the direction of 
surface shear, is denoted 6. This angle is a measure of the 

magnitude of the skewing in the Z' direction of the flow with the 
boundary layer and is a measure of how three-dimensional the 

viscous flow is. All of the viscous calculation methods are uned 


to compute 

the 

growth 

of the 

boundary 

layer, 

which is 

characterized 

by 

four 

parameters . 

These 

am the momentum 

thickness in 

the 

mainflow 

direction, 

T E , the 

boundary 

layer shape 


factor in the mainflow direction, H, the limiting wall streamline 
angle, 6, and the viscous displacement thickness, .>* . 


V 
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Z' 



FIGURE 27. GEOMETRY WITHIN BOUNDARY LAYER 


Along the attachment line the laminar flow equations at 
point i are: 


.407 C #1 /2 
T E t " Roc U E ” 

* Roc u E , 2 

Cl * (Ux e * v y , - U Eyi > 

U E . = (U , 2 ♦ v.V /2 

V. c t 

Si =0 

H, = 2.b4 

£ 

5 , =0. 

L 


123) 

(24) 

(2b) 

( 2t> ) 

(27) 

(28) 
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The shape factor is chosen arbitrarily. The derivative of these 
equations with respect to g\ must be established. As an example, 
consider equation (23): 


3T E^ T E £ 



3 Ue } 

( 29 ) 

3 9k 2C ^ 

^9k 

% 

39k 


1 equations 

(24) 

and 

(25) and 

noting that U x and 


Vy are determined from the surrounding Ui and Vj_, X', and Y', 


equation (29) can be expressed as: 



(30) 


Substituting 


expressions for 

3 ^ 


the derivatives appearing in 


equation (30) results in < A becoming a vector with each element 

3 9k 

representing the rate of change of Tg with respect to change in 
the kth element of the g array. Differentiating equations 
(2o)-(28) results in: 



3H 

3 5* 

39k 

" 39k ~ 

39k 


(31) 


Each of the equations in the viscous flow computation proce- 
dure is differentiated. Most are successive application of the 
chain rule, illustrated above. Therefore, the mathematical 
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expansions for the remaining equations will not be presented, for 
clarity, unless the procedure employed differs substantially from 
the above example. 

Transition along the attachment line is specified by the Rey- 
nolds number based on momentum thickness, Rf, being greater than 
100. Turbulent attachment line flow development is calculated 
from a linear curve fit to the • 0 solution of Smith's method 
(Reference 22): 


.0010925 cj ♦ 83.75 

Te * £ 7 "% 

8 £ ■ 0 


(32) 

(33) 


H. - 1.45 (34) 

5* - 0. (35) 

1 


The value for the shape factor, H., is chosen arbitrarily. 

The attachment line values of the viscous parameters are 
used as starting conditions for the remainder of the boundary 
layer computations . The viscous geometry on the wing is depicted 
in Figure 28 for the example case of fully attached flow, 
attachment line at the leading edge, and transition to turbulence 
along the attachment line halfway across the span. Each of the 
same type of viscous regions across the span is divided into a 
user supplied number of boundary layer calculation points . Fach 
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of these points correspond to the intersection of lines along the 
panel centers and lines of constant Z, . The local nonorthogonal 
coordinate system is based on lines of constant y and constant 
with the angle between the two denoted X. 



FIGURE 28. VISCOUS REGION GEOMETRY FOR A SWEPT, TAPERED WING 

The generalized nonorthogonal coordinate system introduces 
metric coefficients into the viscous equations. These coeffi- 
cients (h^, h2) represent the coordinate stretching factors, 

where an element of length ds on the body surface is given by: 

ds^ = h^ ax + h 2 Ey + h^h^ COS * dx 


V 
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hi and h2 are dependent on X' and Y' and are described in 
Appendix I . 

The flow in the laminar or turbulent reyions is calculated 
by integrating the equations of motion in the X' direction 
(attachment line to trailing edge). The slopes of the T E , B, and 
H terms are calculated at the latest points at which these terms 
are known from equations of the form: 


3 

j£i Ck 3 = (37) 

Cfcj = Coefficients based on the values of Tg, B, 
H, U, V h 1# h2, X' , and Y' . 

, dv , dT E d3 dH 
' dx 1 j = dx 1 ' dx 1 ' dx 1 D= 

RHS^ = Right hand sides, based on values of T E , 

3, H, U, V, hi, h 2 , X' , and Y* . 


Both the laminar and turbulent flow equations are of the form of 
equation (37); however, the laminar flow equations, being a small 
crossflow method, reduce to a 2x2 system of linear equations. 
The solution to equation (37) is established by simple linear 
algebra: 


,dv_ 

'dX 1 


) = 

j 


: cji 

k=l J 


RHSfc 


(38) 


v 
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Trapezoidal rule integration is performed resulting in the 
projected value of T E , 8, and H at point l +1: 


dv 


v i«,i - v i.) + “'i 


( 33 ) 


The displacement thickness does not appear explicitly in the equa- 
tions of motion in three-dimensional boundary layers. Rather, 
the displacement thickness slope is calculated from the 
continuity equation and is of the form: 


dfi 

dX 1- = 


(40) 


where 


dT E 

RHS 4 = Right hand side, based on values of 
d 8 dH 


dX^ dX 


-, 5 *, U, V, hi, h2, X' , and Y' . 


This equation is also integrated in the X* direction 


5 i+i = i 


(41) 


This equation is only applied in regions where the displacement 
thickness is modelled. 


Equations (37) -(41) must be differentiated with respect to 
gfc. Working backwards: 
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H±l . 2S . ,ds* 

3gk 3g* "3x"J ^gk 


d< * 

aiTxn* 


3 5* 

- — - is known from the previous upstream point calculation. 
3 gk 

3( AX', ) 

*■— is available from differentiating equation (20). 

:<gk 

■a ^ r# 

t*— (^£r) is determined from equation (40): 

3gk 3x i 



_i_ (di!) , 

3g k V dX*\ 

, -J- 
3gk 

(RHS 4 ) , 



(43) 

( RHS 4 ) i 

dT E 
f ( dx T ' 

l£_ 

dX* ' 

dH 
X 1 

, 5*. U, 

V, hi, 

h2, X*, Y' ) t 

(44) 

3(RHS ) 

3f 

dT, 

3( — 

h 

3f 
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+ 
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3X ‘ 

3gk 
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3 5 * 

The — s— r- 

3 g* 

calculation 


is available from the previous upstream point 


calculation. y— ^ and are available from the inviscid 

3 h^ 3 h 2 3 y • 

differentiations. T^' and are determined from X , and 

oc * 

is available from differentiating equation (20). The 

derivatives of the slope of T E , £ , and H are determined by 
differentiating equation (3ti): 


5gk ]il Cl 3 ( ^x r) j " 3< 
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* 3 3Cj.j ^ dv 

>1 (Ci 53 9 k l 3X' , j + 39k ‘>X'Y 

dv 3 _ x 3RHSi 3 3Cii 

*ik ( dX r) j “ l-l ° 31 ( "^gk j-1 3 9k 


3RHSi 

39k 




(47) 


(48) 


The Cji‘ 1 *:id ( ^ t ) ^ terms are known from equation (38). The 
derivatives of RHS^ and C^j are determined much the same as 

accomplished for RHS 4 in equation (45). Equation (39) is also 

^ 3T E} 3Hj, 

differentiated to establish expressions for and 

Each of the terms in the mathematical expressions are stored 
as matrices, A, B, C, and D, defined by the following 

expressions : 


N unk 

dU^ - Z dg k 

k -1 

N unk 

dv l * I B tk d ^k 

k -1 

N unk 

d *L ■ I c «k d 9k 

k -1 

N unk 

d 5^ = I D d9k 

k -1 


(49) 


(50) 


(51) 


(52) 


These matrices are retained for use in the viscous-inviscia 
interaction procedure. 


The skin friction at each boundary layer calculation point 
Cf ■ , is calculated from expressions found in the selected 
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methods. The wall shear stress acts on the surface in the 

direction of the limiting wall streamline. Thus, the skin 
friction force acting in the chordwise direction tangential to 
the panel surface is : 

Cfx ■ ■ Cf coe ( Y + 6) (53) 

These values are converted to mean values on each surface panel 
by length weighting. 

The forces and moments on the aircraft surface are deter- 
mined by integrating the pressure coefficient, c p^» over the 
aircraft surface and the skin friction over the surface of the 
wing. A test for convergence based on a .01 change in the wing 
lift coefficient is made. 

The viscous-inviscid interaction is explicitly modelled by 
finding the set of dg^ that most nearly satisfies to first order 
the conditions of the theoretical model. The theoretical model 
can be mathematically stated as follows: 


Aft Attached Flow: 

- Z' 

(54) 

Separated Flow: 

C Pj, * C Psep 

(55) 

Wing/Trailing Edge: 

C Pupr C Plwr 

(5b) 

Fore Wake: Cp 

S i 

* C PTE + S TO t (C PWTE C PTE ) 

(57) 

Aft Wake: 

C Pupr i * C Plwr £ 

(58) 
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where S is surface distance from wing trailing edge, &p 0 T 18 sur- 
face distance from wing trailing edge to fore wake trailing edge, 
Cp TE is average value of Cp at wing trailing edge at each span, 
Cp WTE is value of desired Cp at fore wake trailing edge at each 
span (input by user), and Cp #ep is the Cp at each of the 
separation points. 


A first order mathematical expansion of equations (54-58) 
results in: 


Noting that: 


5*^ + do * - Z ’ 4 ♦ dZ 1 £ 

(59) 

C Pl + dC Pi ” C Psep * dC Psep 

(00) 

C Pupr + dC Pupr C Plwr * dC Plwr 

(01) 

C P t + dC Pl - dC PTE 

(02) 

C_ + dC n - C n . + dC D1 

Pupr^ Pupr i Plwr^ Plwr^ 

(03) 

3Cp -2 , 3U 3V » 

— * - (U V -^— ) 

3 9k 3g)< 3gk 

( o4 


and substituting equations (49)- (52) for the aerodynamic perturba- 
tions, and derivatives of equation (21) for dZ ' , equations (59)- 
(b3) can each be written in the form: 


n UNK 

£ i - I E ik d 9k * 0 
k=l 


(os) 
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where is the difference in desired conditions, e.g., e £ ■ 

Z' - .5*;, and the remainder of the equation is just the first 

i 

order change in e, e.g., 

n UNK 

d ££■ dZ'j - d<5* ■ - i E ^9k (66) 

k-1 

The appropriate equation for each viscous region type is 
applied at each viscous matching point. The solution to this set 
of linear equations would represent the changes in the 

displacement and slope at the region endpoints to satisfy the con- 
ditions of the theoretical model. However, additional equations 
are necessary to establish the X' location of each region end- 
point. These equations, called region endpoint constraints, are 
necessary at each region endpoint that represents an independent 
unknown in the g array. These are the attachment line, 
transition location, and separation location. 

The desired condition at the attachment line is that the 
tangential velocity on the chordwise direction, U, be zero. 
Thus, the perturbation form of the equation is: 

U ATTCH + dUftTTCH * 0. (67) 


The transition constraint equation for eitner the upper or 
Lower surface is based on four possiDilities : fixed by user, 
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4 mainflow transition, crossflow transition, or laminar separation. 

If the user elects to fix transition at a constant local chord 


fraction, 

the 

fraction 

is 

automatically 

converted to an X* 

location, 

X'tr. 

If X'i 

is 

the value of 

X* at the transition 


location initially, the perturbation form of the constraint 
equation is: 

X'i ♦ dX'i - X'tr (*>«) 

Mainflow transition is indicated by the Reynolds number based on 

momentum thickness, R/j, being equal to a curve fit value of criti- 

cal momentum thickness Reynolds number, R^»i>- Rtt i» a function 
of Reynolds number based on surface distance Rx and the 
constraint equation is of the form: 

R>T +■ dRT * Rtt ♦ dRTT (b9) 

Crossflow instability is indicated by the crossflow Reynolds num- 
ber based on boundary layer thickness, Rn , exceeding a critical 
value. Since a change in the crossflow transition point would 
require a change on the geometry in the y direction, this 

equation cannot be applied explicitly. Rather, the most upstream 

location of calculated crossflow transition, if indicated, for 
each span strip replaces X'tr in equation (b8). Laminar 
separation is indicated by the mainflow laminar shape factor, H, 
exceeding 4.0. The constant equation representing the condition 
is : 


H + dH * 4.0 


(70) 
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H + dH - 2.0 (71) 

Equations (68)- (71) are also converted to the form of equa- 
tion (65) by substitution cf aert*3ynaraic perturbations and geo- 
metric perturbations. The number cf equations exceeds the number 
of unknowns. Thus, a method of least square errors is applied to 
obtain the solution g array. 

The solution procedure is to minimize an error function F, 
where: 


F - IWTjJ ej.-2Eijdgj) 2 ♦ l WTfcUk-^kjdgj) (72) 

Z',9' Equations X' Contraints 

The first term in equation (72) represents the appropriate equa- 
tion from equations (59) - (63). The weighting of each equation, 
WT^ , is based on surface distance and user supplied region 
weighting, WTp: 

WT i - (X'i+l - X'i) WT p (73) 

The region weight is used to force the boundary layer type equa- 
tions, equation (59), to be as important as the inviscid 
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4 equations, equations (60)- (63). Normally, a 10,000 to 1 ratio is 
satisfactory . 

The second term in equation (72) represents the reyion end- 
point constraint equations, equations (68)- (71). The weight WT^ 
is calculated automatically in the method and is used to keep f .1 
of these equations more important than equations in the first 
term. 


The method for solving equation (72) is to differentiate F 
with respect to each unknown, dg, and set the resulting terms to 
zero. This converts F to a system of N UN k linear equations, 
which is solved by standard linear alyebra. 

This procedure results in a calculated vector dg*. The y 
array may then be updated: 

9‘k = 9k + d 9k ( ? 4) 

The calculation procedure is then returned to the calculation of 
the A Z ' s . 

This iterative procedure is continued until convergence is 
obtained. This procedure is employed typically over a range of 
angles-of-attack in two degree increments. The last converged 
solution geometry provides the initial guess for the succeeding 
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4 angle-of-attack. Convergence is achieved normally in two to 

three iterations. 


EXAMPLE CALCULATIONS 


To be able to assess the accuracy of the jet-air frame-wing 
viscous effects predictions, several wing alone calculations were 
performed prior to analyzing the YAV-8B powered model. These 
example computations were supported by the MCAIR 1982 Independent 
Research and Development Program. However, these solutions are 
presented here for completeness. The wing alone geometric 
parameters chosen for analysis are shown in Table I along with 
the YAV-8B supercritical wing. Each of these solutions is 

discussed. 


TABLE I. WING GEOMETRIES ANALYZED 


WINQ SECTION 

ASPECT RATIO 

LEADING EDGE 
SWEEP ANGLE 

TAPER RATIO 

TWIST 

(d«g) 

NACA 4412 

6.0 

0 

0 

0 

NACA 0012 

5.5 

20.00 

0 

0 

NACA 64-, -212 

6.0 

37.25 

0.5 

0 

YAV-8B SUPERCRITICAL 
(0.11 t/c ROOT, 0.07 t/c TIP) 

4.0 

36.00 

0.3 

-8.0 


NACA 4412 Rectangular Wing 

The first wing alone geometry analyzed was a rectangular 
planform, NACA 4412 section, aspect ratio o wing. Since 

\ 
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experimental force and moment data were available for this wing 
at midspan, it was possible to assess the accuracy of the method 
without having to deal with strong three dimensional boundary 
layer effects. The predicted force and moment at midspan are com- 
pared with experiment (Reference 27) in Figure 29. The predic- 
tions agree well with experiment. Furthermore, this section 
geometry was analyzed by the Stalled Airfoil Analysis Program 
(Reference 7) with results depicted in Figure 30. Comparison of 
Figure 29 and 30 show that the present method exhibits the same 
trends as the two dimensional method. This result provides confi- 
dence in using the polynomial curve fits to represent the bound- 
ary layer thickness and the Perturbation Analysis Method to 
calculate the inviscid parameters, since the Stalled Airfoil 
Analysis Program does not use either of these approximations . 

NACA 0012 Swept Wing 

The next case analyzed was a 20* swept, NACA 0012 section, 
aspect ratio 5.5 wing on a wall. The predicted forces and mom- 
ents are depicted in Figure 31. While experimental data in the 
form of pressure distributions along lines normal to the leading 
edge are available (Reference 28), overall force and moment data 
are not available. The only force data found in this data set 
was the normal force, C^, on an unswept, NACA 0012 section, 
aspect ratio 6 wing. Thus, this wing geometry was analyzed and 
the predicted normal force is compared with experiment in Figure 
32. Also shown is the predicted normal force for the swept wing 


60 






4 


original 

OF POOR 


PAGr i s 

Quality 


immunnra^ 

fljsfix 10® 


Stalled airtoil analyaia 
program 

O Expanmantai (Mat 27 ) 




- 0.12 



FIGURE 30. FORCE AND MOMENT PREDICTIONS 

NACA 4412 Airtoil 


While analyzing the swept wing case, it was discovered that 
the laminar boundary layer calculations tended to become 
numerically unstable near the suction peak at high angles-of- 
attack (> 15°). It is believed that the primary cause for the 
instabilities was that the transition constraint equations did 
not move the transition location far enough forward to preclude 
the presence of very strong adverse pressure gradients within the 
laminar viscous regions. Successful analysis was accomplished by 
considering the flow to be fully turbulent from the attachment 
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FIGURE 32. COMPARISON OF PREDICTED NORMAL FORCE COEFFICIENT WITH 
EXPERIMENT FOR WINGS WITH NACA 0012 SECTIONS 


NACA b4]_-212 Swept and Tapered Winy 


The last wing alone case analyzed was a winy with an aspect 
ratio of 6, 37.25° leading edge sweep, 0.5 taper ratio, and NACA 
64]_-212 sections normal to the leading edge. The geometry is 
shown in Figure 33. Comparison of predicted forces and pitching 
moment with experiment are depicted in Figure 34. Separation 
occurs at an angle-of-attack of twelve degrees. The preaictions 
do not fully account for the detrimental effects of the flow 
separation. However, comparison of the viscosity included 
predictions with the viscosity ignored predictions indicate that 
a substantial portion of the viscous effects are predicted. 
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SIDE VIEW 

FIGURE 33. GEOMETRY Or SWEPT AND TAPERED WING 


YAV-8B Powered Model 


The final test case, the YAV-8B powered model, was 
accomplished to assess the accuracy of the complete jet-airframe 
interaction method. One of the most important interference 
effects in V/STOL aerodynamics is the lift increment or decrement 
due to power effects in the transition flight regime. If the 
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strong interaction effects are properly modelled, this 
interference can be accurately predicted. 



FIGURE 34. COMPARISON OF PREDICTED FORCE AND MOMENT WITH EXPERIMENT 

FOR A SWEPT AND TAPERED WING 

The YAV-8B powered model with flaps deflected 50* was 
analyzed by the MCAIR jet-airframe interaction method for both 
power off and power on. Viscous corrections were not included in 
these calculations due to the presence of part span flaps, which 
SWAP cannot analyze. The panelled representation of the geometry 
with jets is shown in Figure 35. 




FIGURE 35. PANELED REPRESENTATION OF YAV-8B POWERED MODEL 


The predicted lift coefficient at various angles-of-attack 
is shown in Figure 36 for the YAV-8B with power off. Also shown 
is the range of experimental data. While the slope of the lift 
curve is nearly correct, the level is off significantly for not 
al owing for viscous flow separation on the flap. The ultimate 
goal of SWAP is to calculate this lift decrement; however, 
further developments are required. 
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FIGURE 36. YAV-6B POWER OFF 


Two different approaches were used to predict the lift 
coefficient versus angle-of-attack for the YAV-8B with power on 
(Figure 37). A well designed V/STOL aircraft, such as the 

YAV-8B, uses the large negative angle-of-attack induced by the 

jet-entrainment to keep the flow on the flap attached. Thus, an 
inviscid method should be able to predict the lift curve. The 
first approach used was to utilize a single calculation of the 
jet effects at an injection angle corresponding to U* 
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FIGURE 37. '.AV4B IN TRANSITION 


angle-of-attack. The solid line in Figure 3 7 is the predicted 
result. The predicted lift coefficient at 0 ' angle-of-attack is 
accurate, but the slope of the lift curve is poor. The reason 
for this descrepancy is that as the angle-of-attack is increased 
the jet modelling is not correct. A remedy to this would be to 
recalculate the jet properties at each angle-of-attack. However, 
an alternative to this expensive procedure is to superimpose the 
power off calculated slope of the lift curve starting at the 
predicted lift coefficient at 0* angle-of-attack. This procedure 
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results in the dotted line in Figure 37, which is quite accurate. 
The good results are not surprising since the jets act primarily 
by changing the effective flap deflection, which does not change 
the slope of the lift curve. 

CONCLUSIONS AND RECOMMENDATIONS 

The Subsonic Potential Flow Program is accurate and 
efficient for the analysis of wings, wing-fuselage combinations, 
and airframe- jet configurations. The Adler-Baron Jet-in- 
Crossflow Program accurately predicts jet properties for single 
jets and, when used with the Wooler method for tandem jets, is 
accurate for multiple jets. The Perturbation Analysis Method 
accurately predicts inviscid aerodynamic properties for large 
perturbations in wing geometry characteristic of strong viscous 
effects. The Stalled Wing Analysis Program, incorporating the 
Perturbation Analysis Method and the Wing Design Method, 
correctly accounts for viscous effects for attached flow but nas 
yet to achievs the desired reliabiliry for massive turbulent flow 
separation on wings and plain wing trailing edge devices . 
Efforts are continuing to eliminate this deficiency. 

The present method utilizes the pilot code of the MCAIR 3-D 
panel method, which requires an unnecessarily complex set of 
input parameters. The production version of the MCAIR 3-D panel 
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method has eliminated the tedious input, and, when used with 
available graphics technology, the input procedure is very 
simple. Thus, improvements in the Jet-Airframe-Viscous 
Interaction Method can be expected when the production version of 
the MCAIR 3-D panel method is incorporated. 


During the course of this investigation several deficiencies 
were identified in available analytical methods and experimental 
data. Whereas eliminating these deficiencies represents 
important research topics, it is beyond the scope of the present 
effort. These recommended research goals are summarized below. 


The 'italled Wing Analysis Program uses a small crossflow 
laminar boundary layer method and a 2-D empirical correlation for 
mainflow transition. These methods do not perform well together 
at high angles-of-attacK for swept wings. It is recommended that 
an improved method for 3-D transition be identified or developed 
and incorporated into the method. Futhermore, it is recommended 
that a fully 3-D laminar boundary layer method be tested against 
the existing method to establish the importance of the crossflow 
terms within the laminar region. The fully 3-D turbulent 
boundary layer method used introduces a significant amount of 
numerical processing and complexity. It is not clear that the 
crossflow terms are impor***-*’ enough to warrant the complexity. 
Thus, it is recommended that a small crossflow turbulent boundary 
layer calculation routine be incorporated into the method and 
compared with the fully 3-D method. 
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A search of the literature for experimental data for 
comparison with prediction indicated some short comings. First, 
very little pressure data on swept and tapered wings character- 
istic of modern aircraft were found. Second, pressure data on 

wings with flaps were not identified. Third, pressure data on 

wing-fuselage combinations with and without strong jet 
interactions were scarce. Lastly, detailed pressure data in the 
wake of jets for realistic configurations were not found. It is 
recommended that detailed experiments be performed to enlarge the 
data base to include the identified cases. 
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APPENDIX I - METRIC COEFFICIENTS 

The metric coefficients (h^, h 2 » h 3 ) establish the relation- 
ship between the Cartesian coordinate system (X, Y, Z) and the 
nonothogonal boundary layer coordinate system (X‘, Y' ) on the 
surface of the wing. The arc length on the surface of the wing, 
ds, is determined from: 

ds 2 = h 2 ! dx 2 + h 2 2 dy 2 + h 3 2 dxdy (1-1) 

where 

h 3 2 - 2h]_ h2 cos X (1-2) 

and X is the angle between the X' and Y‘ axes. These definitions 
arise out of the use of the metric tensor, which also establishes 
the relationships between the two coordinate systems: 


h l 

h2 

h3 
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( T3c r) 
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3 X ' 
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(1-3) 

(1-4) 

(1-5) 


Thus, if analytical 
(X', Y'), equations 
and h 3 . 


relationships are known between (X,Y,Z) and 
(I-3)-(I-5) can be used to calculate h^, h2 
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Consider the wing as broken into two surfaces at the attach- 
ment line, and one of these surfaces has been unwrapped and laid 
flat in a plane. Assume that the unwrapping is accomplished in a 
way that keeps the lines in the X 1 direction continuous at the 
leading edge. 

The (X', Y' ) origin is located at the junction between the 
win:; leading edge and wing root. The functions of Y', a, b, c, 
and d, are defined as (see Figure Al): 


a(y' ) * X 1 tt.f 

(y 1 ) 

(1-6) 

b(y ' ) = X' ATT 

(y ' ) - x * ile (y * ) 

(1-7) 

c(y ' ) - X' TR 

(y' ) - X 1 !££ (y* ) 

(1-8) 

d(y ' ) * X'te 

(y ' ) - X'-t-r (y ' ) 

(I-J) 


From these definitions and by considering the wing to be flat, X, 
Y, and Z can be determined as: 


Laminar Region: X = a(y') + b(y') + ( a (y ‘ = U" ) ~b(y " ‘ -u ) ^ (x'-b(y'-o)) ( 
Turbulent Region: X * a(y' ) + c(y') + d ^y¥,Q r (x'-c(y’»U)) (I -11) 


Y * 2 

(1 + a* 2 ) 

Z = 0 

where a' is the derivative of a with respect to y' 

, _ da (y 1 ) 
dy' 


(1-12) 

(1-13) 

(1-14) 


1 - 10 ) 
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VIEW AA 


FIGURE A1. DEFINITION OF GEOMETRIC PARAMETERS THAT ESTABLISH 

METRIC COEFFICIENTS 


This definition of Y may seein to be unnecessarily compli- 
cated ’ :ause of the obvious trigonometric relationship between Y 
and Y' . However, the viscous computation methods require the 
local value of X and Y' at each calculation point. Since the 
region boundaries at each span station are allowed to move 
independently, these local values can vary from one point to the 
next . 

The angle between the local coordinates, Xj_, is defined as: 


tan 



(I-1S) 
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Thus : 

Laminar Region: tan 


»■ + h« + l a 1 - b’ W y 1 -h ( y 1 -G n 
«(y'"0)-b(y '*0) 


(I-lb) 


Turbulent Region: tan X^ 


♦ c' + _11_ (x'-cty'-O)) 

d ( y 1 *0 ) 


(1-17) 


The metric coefficients hj_ f h2 and h 3 are calculated by differen- 
tiating equations ( I —10 ) — (1-13 ) and substituting into equations 


h3 


The 

results for 

hi * 

d(y') 

d(y'-O) 

h2 * 

'[(a 1 + c* + 


d(y') a' + c 1 


(1-18) 


d' 


d(y ’-0) (x‘ -c(y'»U) ) ] + 1 


1 + a’ 2 


(1-19) 


d(y'-O) 


(X + a' 2 ) 1 ^ 2 


The metric coefficients for the laminar regions are not 
used. The gradients of h]_, h 2 » and h 3 in the X' and Y' 
directions are established by differentiating equations 

(I-l8)-(I-20) . Each of the equations presented in this appendix 
are differentiated with respect to the g^ array for inclusion in 
the viscous inviscid interaction procedure. 
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